Interface Motion and Pinning in Small World Networks 
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We show that the nonequilibrium dynamics of systems with many interacting elements located on 
a small-world network can be much slower than on regular networks. As an example, we study 
the phase ordering dynamics of the Ising model on a Watts-Strogatz network, after a quench in 
the ferromagnetic phase at zero temperature. In one and two dimensions, small-world features 
produce dynamically frozen configurations, disordered at large length scales, analogous of random 
field models. This picture differs from the common knowledge (supported by equilibrium results) 
that ferromagnetic short-cuts connections favor order and uniformity. We briefly discuss some 
implications of these results regarding the dynamics of social changes. 



Small-world networks have received a great deal of at- 
tention in the past few years, in particular for their real- 
istic description of the topology of the interactions that 
take place among populations in various biological, so- 
cial or economical systems [1]. An important feature of 
small worlds, that is not shared neither by regular lat- 
tices nor random networks, is the interplay that exists 
between local (or "physical") interactions, e.g. between 
nearest- neighbors, and non-local ones, involving nodes 
(or agents) separated by large distances but connected 
through short-cuts. Among other outstanding topologi- 
cal properties, the effective space dimension of such net- 
works grows linearly with their size [2] , even if the frac- 
tion of sites with short cuts is very small. 

The strong connectivity of small worlds usually en- 
hances dramatically cooperative effects, as predicted by 
epidemic models of spreading of diseases [4] , or of prop- 
agation of conventions or rumors in social systems [5]. 
Naturally, many models of social dynamics have been 
inspired from the Ising model [6]. The Ising model on 
a small world exhibits ferromagnetic order at low tem- 
peratures even in one dimension [7], while, in higher di- 
mensions, the critical temperature is increased compared 
with that of the regular lattice [8] . In addition, the fact 
that the transition is of mean-field nature agrees with the 
intuitive argument that each site is effectively close to a 
large number of sites due to the short-cuts of the lattice. 

However, because of their inherent random topology, 
one may ask whether in some situations small- world net- 
works would not rather exhibit features characteristic of 
disordered systems. In this work, we study as a basic ex- 
ample the nonequilibrium dynamics of the Ising model, 
as observed after a rapid quench from the high tempera- 
ture phase to the ferromagnetic phase. We show that the 
random (all ferromagnetic) connections that enhance or- 
dered states at thermodynamic equilibrium, are respon- 
sible in the present case for very slow dynamics and sta- 
bilize at large times configurations that, instead of being 
uniform, are spatially heterogeneous. At zero tempera- 
ture, systems do not perform long range order dynami- 
cally, but remain asymptotically trapped in metastable 
states characterized by a finite domain size. These fea- 



tures are reminiscent of nonequilibrium processes in the 
random field Ising model (RFIM) [9] , in binary mixtures 
with fixed impurities [10], as well as in a few social mod- 
els on regular lattices like the voter model [11]. This has 
to be contrasted with the much more efficient phase or- 
dering kinetics of the Ising model on regular lattices (or 
Model A in the lexicon of Hohenberg and Halperin [12]), 
where the mean size of ordered domains grows with time 
as i 1 / 2 [13]. Our present analysis focuses on the motion 
of domain walls between "up" and "down" domains, and 
shows evidence of competing effects between surface ten- 
sion and pinning (or localizing) effects. 

We use a standard model of small-world network [14] 
consisting of a regular square lattice (or a chain in ID) 
composed of N nodes connected to their nearest neigh- 
bors. For each site, we then establish with a probability p 
an additional connection, or short-cut, linking the consid- 
ered site to an other site chosen at random in the lattice. 
(We do not remove the nearest neighbors connections.) 
For p — 0, the lattice is regular, while for p = 1, the net- 
work is strongly disordered. Here, we will consider only 
the so-called "small world" limit, that corresponds to the 
case p <C 1 , where connections are mainly local and only 
long-ranged for a small fraction of nodes. 

On a fixed network, we then assign to each node a 
spin- like variable Si = ±1: it represents a social con- 
vention, initially chosen at random for each node. At 
each time step, each node updates its convention in or- 
der to reach a better consensus with the nodes it is con- 
nected to. In other words, the system follows a zero 
temperature Glauber dynamics with the Hamiltonian 
H = JijSiSj, where the sum is performed over 

all possible pairs of nodes. = 1 if sites i and j are 
connected, = otherwise. At each step, a spin is thus 
chosen at random and flipped. If H decreases, does not 
change, or increases, the new configuration is accepted 
with probability 1, 1/2 and 0, respectively. 

In regular networks (p = 0), the system evolves to- 
ward a minimum of H (all S^s equal to +1 or — 1). 
Transient configurations are characterized by the pres- 
ence of growing and competing ordered domains of "up" 
and "down" spins. The large time dynamics is controlled 
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by the motion and annihilation of interfaces (or domain 
walls) that separate these domains. As for many other 
systems ordering in phases with broken symmetries, its 
time evolution is self-similar: The two-point correlation 
function, C(r,t) = (Si(t)Si+ r (t)), obeys a scaling rela- 
tion C(r, t) = f(r/£(t)), where / is a scaling function, 
while £, the domain size, grows as t 1 / 2 [13]. 
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FIG. 1. Asymptotic correlation length £o in lattice spacing 
units as a function of the reconnection probability p, in ID 
and 2D. I? j L x is the interface characteristic length in 2D. 

On the contrary, in small- world networks (p ^ 0, 
p <C 1) one observes after some time that the typical 
domain size of ordered spins saturates to a finite value, 
which decreases when the density of short-cuts (p) in- 
creases. For a one dimensional chain (of length L = 10 5 ), 
we plot £(t = oo) as a function of p in Figure 1. The cor- 
relation length is determined from the half-width of C 
averaged over 10 networks and initial conditions. A be- 
havior £(t = oo) oc 1/p can be observed, l/p represents 
the characteristic size of the one dimensional network, 
i.e. the average distance between two nodes that have 
long range connections (or "influent" nodes). 

Influent sites strongly affect the motion of interfaces. 
At low p, most of these nodes are characterized by one ad- 
ditional connection. On Figure 2a, two nodes far apart, 
A and B, are connected, and Sa = —Sb- Any inter- 
face I passing through node A leftward can not jump 
back toward the right, since it is energetically unfavor- 
able. Therefore, at large times, through interface motion, 
influent nodes will tend to be (irreversibly) connected to 
nodes that have the same spin (a situation analogous to 
assortative mixing, as observed in some real life networks 
[3])- 

This argument can be extended to a succession of do- 
mains. Figure 26 illustrates a typical large time configu- 
ration: In this example, interface I\ stands between two 
influent nodes A and B with opposite spins: for the rea- 
son mentioned above, I\ can not jump to the left of A, 
nor to the right of B. The interface is then localized, i.e. 



restricted to perform a random walk within the interval 
L4,£?]. Interface I\ is therefore unable to annihilate with 
interface I2, that is localized between C and D. Hence, 
the two disjoint black domains can not merge to form 
a bigger one, what would happen in the standard Ising 
model. The domain size, or correlation length £, does not 
exceed the distance, of order 1 /p, that separates influent 
unlike (antagonist) successive nodes. 
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FIG. 2. Up and down domains. Domain walls I n become 
localized in ID. 

We also find that the structure of frozen configurations 
obeys a scaling relation with p, i. e. that they are statis- 
tically independent of p via proper length rcscaling. The 
asymptotic correlation function C(r,t — 00) is plotted in 
Figure 3 as a function of the reduced variable r/^ 00 (p), 
for various values of p. Data collapse rather well on a 
single curve. At short times, the kinetics is not affected 
by the small world structure of the lattice, and that 
starts growing as t x l 2 . When interfaces become local- 
ized, the structure can be roughly seen as the one given 
by the standard Glauber dynamics of the ID Ising model 
stopped at a time (1/p) 2 . Since that problem obeys dy- 
namical scaling, frozen configuration should scale with 
parameter p. However, this picture is not quantitatively 
correct, as the scaling functions in both problems slightly 
differ. 

Finite, low temperature effects are quite subtle in one 
dimensional small-worlds since they do not destroy the 
ferromagnetic order observed at T = 0, unlike in usual 
Ising chains. One can interpret here the order/disorder 
transition temperature T c as the temperature where or- 
dering via interface jumps over localizing barriers (that 
enable further domain merging) no longer overcomes 
disordering happening within domains (subdomain cre- 
ation). In Figure 26, interface I2 (or I{) can jump in 
interval [B,C] at a rate r a = exp(— 2/T). Besides, the 
rate at which any spin among the p~ x spins of interval 
[D,E] would flip is r& = p^ 1 cxp(— 4/T): it is roughly the 
rate at which a white domain is created and can start to 
grow. Qualitatively, the order/disorder transition occurs 
when r a = rv This gives: T c ~ —2/ lnp, a expression de- 
rived (with a oc sign) in ref. [7] using the replica method. 
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The above relation may be exact as p — > (as the numer- 
ical prefactors in the different rates become irrelevant). 
Simulation results (not shown) give T c = — 2.3/ In p for 
p = 0.03. 



p = 0.2; 0.1; 0.07; 0.05; 0.03; 
0.02:0.01:0.005 
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FIG. 3. Two point correlation function in ID (dotted lines) 
and 2D (solid lines) as a function of r/^, for various p. 




FIG. 4. Frozen domains (in gray) for p = 0.05 (6% of the 
system total area). The white and black dots represent the 
"like" and "unlike" influent nodes, respectively. 

In two dimensions, one also observes that random ini- 
tial configurations freeze at large times. Figure 1 displays 
as a function of p the asymptotic correlation length (de- 
termined from C(r,t — oo) averaged over 8 networks with 
1500 2 spins), as well as the length associated with inter- 
face density. The latter is defined as L 2 /£, where L is 
the system linear extend and £ the total length of all 
boundaries. Both length scales remain proportional to 
each other when varying p, suggesting that frozen config- 
urations can be characterized by one characteristic length 
scale, referred to as the "domain size" , Roo (p)- Numerical 
results suggest that varies as an inverse power-law of 
p, with an non-trivial exponent close to —2/3 over nearly 
two decades. Surprisingly, R x does not scale as p^ 1 ^ 2 , 
the characteristic length scale of the network [2]. Once 
again, the spin-spin correlation function at t = oo scales 
rather well as C(r) — /(r/£oo(p))> see Figure 3. 

Figure 4 shows a typical frozen pattern at p = 0.05. 
The positions of the "influent" spins are marked by dots: 



in white, those which are connected to an other influent 
spin of same sign ("like" pairs, of number density n{), 
in black those connected to a spin of opposite sign ( "un- 
like" pairs, number density n u ). Initially, ri[ ~ n u , but 
as coarsening proceeds, "unlike" dots turn more easily to 
"like" than the contrary, as in one dimension (Fig. 2a). 
Once again, mixing tends to be assortative (n/ > n u ), but 
with the increase of ni, at some point, there are no more 
possible moves toward better consensus. We find nu- 
merically that coarsening stops and interfaces get pinned 
when ni ~ 1.86 n u . 

The finite domain size can be interpreted as the result 
of competing effects between surface tension (the driving 
force for domain growth) and energy barriers created by 
the multiplication of influent "like" sites. We picture the 
system as a collection of L 2 / R 2 domains of radius R, and 
estimate its energy change when domains coarsen from 
R to R + dR (dR > 0). The usual contribution from 
surface tension is SEi oc —2L 2 dR/R 2 . Meanwhile, the 
number of influent nodes that flip spin is proportional to 
2pRdR(L 2 /R 2 ). "Like" nodes turn to "unlike" (with an 
energy cost per spin of 2), and reversely (with an energy 
decrease of -2). The total energy difference thus reads 
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The second term is positive and dominate at large R. 
Hence, coarsening is arrested when SE = 0, or R x ~ 
p^ 1 . This argument is somehow similar to the (equilib- 
rium) Imry-Ma argument for the RFIM [15]. Yet, an 
important difference is that here the average magnetic 
field felt on influent nodes (or "impurities") is not zero, 
but has been biased (n; ^ n u ) due to previous spin flips. 

The above continuous Imry-Ma-like argument quali- 
tatively explains frozen states, but over-estimates R^ 
(~ p^ 1 instead of p~ 2 / 3 ). The exponent —2/3 can be 
explained as an effect of the square lattice. As shown 
on Figure 5a, a single influent "like" node located at a 
domain corner can disappear through the diffusive mo- 
tion of a step. Figure 56 represents then the simplest 
distribution of "like" nodes such that the hatched do- 
main can not shrink. It is composed of two right-angle 
corners {^,^1,^2} and {B, B\, B 2 } defining a square 
rxr. If the other white nodes {Di, ..D n } contained in the 
square do not form any right-angle corners, then this re- 
gion encloses the smallest (or "critical" ) pinned domain: 
any bubble of hatched region comprised in the square 
and that does not contain both corners {A, A\, A2} and 
{B, Bi, B2] will shrink. Any larger bubble will not. We 
now calculate the probability PfroczcM that a configu- 
ration such as represented in Fig. 56 has a size r, and 
then identify r* such that -Pfreeze^*) is maximal with the 
asymptotic domain size R^ in the disordered medium. 

Given the node A located at the origin, the probability 
that there is at least one white dot (Ai) on the same 
line within a distance r is Pi(r) = 1 — (1 — pi) r , with 
pi jp = nij (n; + n u ) the fraction of influent nodes that are 
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"like" (here, the numerical value of this ratio — close to 
0.65 — is unimportant and could be set to 1). Therefore, 
-PfrcczcM = Pi[Pi(r)] 4 P 2 (r), with P 2 the probability that 
the -D„'s do not form right-angle corners, i.e. that each 
node Di is at least located on a line or a column not 
occupied by an other Dj (see the dotted lines in Fig. 
5b). P 2 can be approximated as 

(r-2) 2 

P 2 (r)^ J2 a-Pi) (r - 2f - n P?Cl_ 2) 4l-Pi(r)] n , (2) 

n=0 

or P 2 (r) — (1 — PiPi{r))^ r ~ 2 ^ ■ In the sum (2) we have 
multiplied the probability of having n white dots inside 
the square by the probability [1 — P\(r)] n that n inde- 
pendent dots have no neighbors on the same line within 
r. For n small, no or few corners can be formed any- 
way, so that relation (2) slightly under-estimated P 2 , 
since a small fraction of empty sites are counted twice 
([1 — Pi(r)] n < 1). For n large, on the contrary, rela- 
tion (2) over-estimates P 2 , since it is impossible to locate 
many dots without forming corners (while [1 — P\(r)\ n is 
small but ^ 0). We suppose that both errors compen- 
sate. This factorization enable us to compute the most 
probable square size r* analytically. P\ (r) increases with 
r, P 2 (r) decreases with r, and Pf Teeze (r) has one single 
maximum. Assuming p « 1, r > 1, rp < 1, we find 

that r* = P[ 2 ^ 3 cx p~ 2 ^ 3 , in agreement with the numeri- 
cal results. 
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FIG. 5. a) Free and b) pinned domain in presence of "like" 
influent nodes. 

To summarize, we have shown on an example that as- 
sortative mixing in small world networks can dynami- 
cally generate frozen metastable states : At large times, 
some influent nodes have simply no immediate interest 
to evolve. These results suggest that long term dynam- 
ics in highly connected social systems can produce spatial 



heterogeneities (or segregation), despite that these con- 
figurations are not the most desired ones by individual 
agents. A similar picture, in agreement with some em- 
pirical observations, was drawn recently from antifcrro- 
magnetic models on scale- free networks [16]. Right after 
strong political changes (in Eastern European countries 
in 1989, in Mexico in 2000) the evolution of reforms can 
be fast, but rapidly social inertia takes over and renders 
further adjustments difficult or null. Physically speak- 
ing, the response of social systems to external forcings 
(i.e. large-scale policies) is susceptible to exhibit some of 
the interesting features known for disordered systems [9] . 

While revising the manuscript, we became aware of a 
similar study on the voter model on small- worlds [17]. 
We acknowledge fruitful discussions with G. Cocho, J. 
Vinals and R. Boyer. 
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